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£>. ■ Abstract 

in 
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For Dirac equation, operator-invariants containing explicit time-dependence in parallel 
(N ' 

r-. to known time-dependent invariants of nonrelativistic Schrodinger equation are introduced 



and discussed. As an example, a free Dirac particle is considered and new invariants are 
constructed for it. The integral of motion, which is initial Newton-Wigner position operator, 



& . 

is obtained explicitly for a free Dirac particle. For such particle with kick modeled by 



delta-function of time, the time-depending integral, which has physical meaning of initial 
momentum, is found. 

1 Introduction 

The usual time-independent integrals of motion of stationary quantum system like energy or 
angular momentum are related to symmetry of the Hamiltonian, which is connected with time 
translational symmetry or rotational symmetry of the potential |IJ. For nonrelativistic quantum 
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systems with nonstationary Hamiltonians (for example, for the parametric oscillator), quadratic 
in position and momentum integrals of motion of classical problem, known long time ago [§], 
were shown to exist for quantum problem as well [|J. A connection of finding wave function of 
nonstationary Hamiltonians with existence of time- dependent integrals of motion has been also 
elaborated [Q. 

For the same problem of nonrelativistic parametric oscillator, linear in position and momentum 
time-dependent integrals of motion have been found [5-7] and have been related to the Green 
function of Schrodinger equation ||. A connection of the integrals of motion with quantum 
propagator turns out to be general property related to Schwinger action principle ||. Futhermore, 
an extension of the Nother theorem from time-independent integrals of motion to time-dependent 
integrals of motion has been done [10-13]. The method of finding time-dependent invariants was 
developed in recent works [14-17]. 

The time-independent integrals of motion as energy and angular momentum are used as well 
in relativistic problems of scalar and spinor particles in the framework of Klein-Gordon or Dirac 
equations []]]. But up to our knowledge till now the time- dependent integrals of motion for 
relativistic quantum systems have not been discussed and there is still no one example of time- 
dependent integral of motion for relativistic quantum system. 

The aim of this work is to find new time-dependent integrals of motion for free Dirac particle 
and to extend the results known for nonrelativistic propagator of Schrodinger equation, which 
may be obtained as eigenfunction of time-dependent invariants || to the case of Dirac equation. 

The problem of position operator in relatvistic domain is more complicated than for Schrodinger 
wave mechanics. There are known results for Dirac equation to introduce an operator playing role 



of position operator |I5 , pTQf| . Newton- Wigner position operator is intensively discussed, in recent 



publications as well [20] . Thus, another goal of our work is to construct a time-dependent inte- 



gral of motion, which coincides with Newton-Wigner position operator for initial time moment. 



Production of electron-positron pairs in heavy ion collisions [21] was studied recently in frame 



of the model of short pulses described by delta-kickes [g^]. Dirac equation was used in which a 
potential was switched in for very short time corresponding to ultrafast interaction of high energy 
of colliding ions. In this connection, it is worthy to investigate what type of integrals of motion 
does exist in such potential. We construct explicitly expressions for some time-dependent integrals 
of motion for Dirac equation with short kickes modeled by delta-function of time. 

The paper is organized as follows. In Section 2, we review the known time-dependent invariants 
for free nonrelativistic scalar Schrodinger particle (readers having interest to detail review we 



address to Refs. [23[ 0] ) and generalize the nonrelativistic invariants to the case of free Dirac 
particle. In Section 3, we discuss the time-dependent invariants and their properties for Dirac 
particle in electromagnetic field. In Section 4, we discuss Newton-Wigner position operator and 
its time-dependent analog, which is an integral of motion for Dirac free particle. In Section 5, 
an extension of the integrals of motion found for free Dirac particle is done for the kicked Dirac 
particle in frame of the delta-pulse model. 

2 Invariants for Free Dirac Particle 

Let us consider a free Dirac equation in standard form 

dip 
i — = (ap + pm)ip ; h = l = c, (1) 



where pi = —i d/dxi (i — 1, 2, 3) are the components of 3-vector p, m is constant mass, and the 



4x4-matrices «j and /3 have the form 



o 



a 
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(2) 



where the vector-matrix a has three components on. The Pauli 2x2-matrices which are compo- 



nents of vector-matrix o are chosen to be 
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v 1 °j 
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03 



-1 



(3) 



The aim of this section is to discuss for a free Dirac particle the existence of a time-dependent 
integral of the motion. 

Let us remind that for a free Schrodinger particle with the Hamiltonian 



H 



2m ' 



Pi = -i 



_d_ 
dxi 



z = l,2,3, 



(4) 



there exist the following invariants: momentum operator described by 3-vector 



Po(t) =p, 



(5) 



which commutes with the Hamiltonian (^), and the initial coordinate vector-operator Xo(t) with 

the components 

/ \ it d . . 

x 0i {t)=Xi^ — - . (6) 

m dxi 

Here Xi are components of 3-vector x. The vector-operators x (t) and po(t) have the property that 



wit) i x (t) i m) = (v»(o) i * i v(o)) 



(7) 



and 



m) i p (t) i m) = (v>(o) i p i ^(o)> 



(8) 



for any solution to the Schrodinger equation. At the same time, these operators may be expressed 
in terms of the evolution operator U(t) of the free Schrodinger particle. This evolution operator 
satisfies the equation 

i% = 2-0 (9) 

dt 2m K J 

with the initial conditions 

U(0) = E. 

The operator E is the identity operator. In terms of the operator U(t), the invariant Xo(t) (§) is 
expressed as follows 

x (t) = U(t)xU- 1 (t), (10) 

and the operator po(t) @ is expressed in the similar manner 

p (t) = U(t)pU- 1 (t). (11) 

The operator U(t) is the solution to Eq. @ and has the form 

U(t) = exp (-it ^j . (12) 

Thus, the exact forms of the operators xq and po are the result of commuting of the operator U 
( |T2| ) either with operator p or with operator x in formulas (|TT| ) and ([T0|), correspondingly. At the 
same time, these operators satisfy the system of equations 

; (is) 

(14) 

For this simple problem, we can use a trick to guess the form of the operators (0) and @ on the 
basis of classical intuition. 



dx 

+ i[H 1 x G \ -- 


= 


dpo ,.,„., 


= 



Let us write down the classical solution for the free motion equation. We have 

p(t) = Po; 

x(t) = x + v t; (15) 

P = 171V . 

Here po and vo(po) are the initial values of the coordinate and velocity (momentum) of the particle. 
Expressing these initial values in terms of the instant values pit) and x (£), in view of (JIB]) , we 
have 

Po = Pit); (16) 

Xo = x -t^-. (17) 

m 

"Quantizing" these initial values we obtain the quantum invariants (^) and (^). 

From formulas ( |TTiD and (|H|), it follows that the propagator of a free Schrodinger particle 
G (x, x', t) satisfies the system of equations 

x {t)G(x,x',t) = x'G{x,x',t); (18) 

d 
p (t)G(x,x',t) = i—G(x,x',t). (19) 

ox' 

Thus, knowing invariants (|5|) and (0) we have the equations for matrix elements of evolution 
operator in the coordinate representation. 

Now let us discuss what properties of the Schrodinger invariants po and xo may be generalized 
for the Dirac particle. It is the general problem. But first we have to consider the problem applying 
it to the simplest case of a free particle. Equation (|l|) is the equation of the same evolutionary 
type as the Schrodinger equation. Thus, the 4x4-matrix Un{t) may be introduced as evolution 



operator. This operator-matrix satisfies the equation 

i °^- = H D U D , (20) 

where the Hamiltonian Hjj is the 4 x4-operator- matrix 

H D = ap + (3m. (21) 

The solution to (|20"D has the same form that in the Schrodinger case 

U D (t) = exp (-iH D t) . (22) 

Let us now generalize formulas ([It]) and ([II]) of the Schrodinger case to the Dirac case. We have 
the operators 

PoD(t) = exp(-iH D t)p exp(iH D t); (23) 

Xooif) = exp(—iH D t)xexp(iH D t). (24) 

Equation ( p3[ ) immediately gives the equality for p D (t) , namely, 

Po D (t)=p, (25) 

as well as it was in the Schrodinger case, since the Hamiltonian (|2~ID of the Dirac free particle 
commutes with the momentum p. In order to calculate expression (|24]) , we could use the Campbell- 
Hausdorf formula 

e B Ae- B = A + [B, A] + i [B, [B, A}} + ■ ■ ■ + i [B, [B[. . . [B, A]]] ...] + •••. (26) 

But since the first commutator (if one takes B = —itHo and A = X{ which is the i-th component 
of the vector x ) 

[B, A] = [-it (a k p k + (3m), x»] = -it {-i5 ik )a k = -ton (27) 



is not c-number, we need to calculate the next commutators in ( |26| ) 



[B, [B, A]] = [-it (a k p k + (3m), -taj] = it 2 ( [a k ai\p k + m[(3ai 



m 



The matrices ckj and (3 obey the known algebra 



a.ia k + a k a.i 



2 Si. 



a>il3 + [3a>i 



0; 



(29) 



«? = # 



1. 



From this algebra, taking into account (EBT) it is not obvious how to use these formulas for calcu- 
lating each term in series (EBT) as well as the total sum. 



To evaluate this sum, we use the specific property of the Dirac Hamiltonian for free particle. 
Let us introduce —it = t. Then the matrix Hp has the form 



H 



( - \ 

m ap 



D 



\ ap —m I 
The square of this matrix is the following diagonal matrix 



Hi 



m 2 + p 2 



e 2 E 



where 



I vn l + p l j 



2 -2 , 2 

e = p + m 



(30) 



(31) 



(32) 



and E is the 4x4-identity matrix. Then one has 



H D = e H D , H n = e E 



, . . . , ±± D 



H 



2n-l J2n-\ n D 



e 



iHl 



, H% = e 2n E. 



(33) 



Thus, 



e 2 r 2 £ 3 r 3 jj fj 

e TH ° = E + teH d + —-E + — + ■ • • = £ cosh ex + — - smher . (34) 



This operator in the matrix form is 



e rH D 



cosh.ET + E 1 msinh.ST e x a/3 sinh £:r 



(35) 



£ -1 0psinh£T cosher — e~ 1 msinheT 

The matrix e~ tHd differs only by sign before the term with sinus hyperbolic. At r = , this 
matrix coincides with the identity matrix E . 

The Green function Gik(x,x',t) is the 4x4-matrix (i, k — 1, 2, 3, 4) containing matrix ele- 
ments of the evolution operator Unit) ■ It is connected by the relation similar to ([T^) and (|T^) 
with constructed invariants Xou(t) and PoD{t) = P which are 4x4- matrices as well. So 

d 
PoD(t) G ik (x, x , t) =i — G ik (x, x , t) (36) 

and 

(x 0D (t) ) iV Gii k (x, x', t) = x'iGik(x, x, t) . (37) 

Here we have in mind summation over two equal indices i' from 1 to 4 and no summation over 
index % . The operators in the matrices xqd and pqd act on the first argument x of the propagator 

^ikyX) X i Z). 

Now we will write down four matrix elements of the matrix xoo{t)- This matrix consists of four 
2x2-blocks. We will label these blocks as (xo)n, (^0)12, (^0)21, an d (^0)22 • ( in the case of the 
free Dirac particle, the matrix structure of the operator Pooit) is simple because it is proportianal 
to unity 4x4-matrix.) Thus, the 2x2-block (xo(t) )n of invariant xqd is the following one 

771 1 

(xo(t))n — cosh et x cosher sinh et x - sinh et 

E E 



+ m 
op 



1 . , . 1 

- sinh et x cosh et — cosh et x — sinh et 

E E 



sinh et x — sinh er . 

£ E 



(3? 



Here 



£ = A/m 2 +p 2 



• 9 r 



1,2,3); T = -it. 



(39) 



We must take into account the order of the operators in the product in this formula because they 
do not commute. The next block (x (t) )i 2 has the form 



(x (t) 



12 



in 



— sinh et x - sinh et sinh et x — sinh et 



+ — sinh et x cosh et — cosh et x — sinh et . 



(40) 



Of course, we can use the commutator [xi, Pk] = iSik , (i, k — 1, 2, 3) in order to rewrite the terms 
in these formulas. For the block (x ) 2 i , we have the formula 



{x (t) 



'21 



III 



- sinh et x — sinh et 

e E 



sinh et x - sinh et 



+ — sinh et x cosh et — cosh et x — sinh et . 



(41) 



The block (0:0)22 of the 4x4-matrix-invariant x (t) has the form 



(^0 (t) ) 22 — cosh et x cosh et — sinh et x sinh et 

E 2 



+ m 



- cosh et x sinh et sinh et x cosh et 

e E 



— sinh et x 



(42) 



From the structure of these blocks, it is evident that it is not easy to obtain the invariant-operator 
x using the Campbell-Hausdorf formula. Thus, we have generalized the integrals of the motion 



for the free Schrodinger particle to the case of the free Dirac particle. By constructing the 4x4- 
matrix-operator p £, and Xqd obey the Weyl-Heisenberg algebra with commutation relations 

[poDk, Xqdi] = -iSik • (43) 

From these operators, it is possible to construct some invariants with different commutation 
relations. So taking the operator 

Li(t) = eu k xo D ipoDk, i, k, I = 1, 2, 3, (44) 

one can check that the three time- dependent operators obey the commutation relations of the 
angular momentum 

[Li(t),L k (t)]=ie M Li(t). (45) 

On the other hand, at t — these operators coincide with the orbital angular momentum operator 
L(t = 0) = [xxp], which does not conserve, but the operator ([44] ) is invariant. Along with the 
operators Pod, the invariants L(t) form the Lie algebra of the generators of the motion group in 
three-dimensional space. For the free Dirac particle, the conserved operator (which is the sum of 
orbital and spin momenta) j = [xxp] + er/2 . It commutes with the Dirac energy operator Hp . 
Due to this, one can write down the equality 

U D L(t = 0)U D 1 + ^U D aU D 1 = L + ^. (46) 

Thus, the orbital angular momentum described by the operator L(t) may be expressed as 

L{t) = L{t = Q) + a -- l -U D aU D 1 . (47) 

The matrix Unit) [c/2] U^ l {t) = S(t) is the following integral of the motion which describes the 
initial spin operator. It may be obtained in the form of 4x4-matrix operator using formulas (f24j) 



and (|38|)-(f4"2D by the simple substitution: x —> a in all these formulas, all other terms in the 
formulas being unchanged. 

3 General Consideration of Operator-Invariants for Dirac 
Particle 

The analysis of the previous section was based on the fact that the free Dirac particle motion 
may be described by the equation of an evolutianary type with the concrete Hamiltonians (f20|)- 
(pnp. After that, we used the free propagator or evolution operator for free Dirac particles in 
order to construct the invariants PoD(t), xon(i), and L(t) which contained dependence on time. 
But it is clear that the same procedure may be fulfilled not only for the free Dirac particle. Any 
evolutianary type equation for bispinor ip of the kind 

t^t = H D iP (48) 



+ {3m + V{x, t) , (49) 



with the Hamiltonian 

Hd = \ a p A(x, t) 

I L c 

in which we have taken into account also the possibility of presence of electric and magnetic fields, 
posesses the evolution operator or propagator Unit), which makes the bispinor %p(t) from the 
initial bispinor ip{Q) satisfying Eq. ( |48|) 



iJ>{t) = U D (t)iJ)(0). (50) 

If for the more complicated cases with electric and magnetic fields the evolution operator TJjj{t) 



is found, we can construct invariants ii)(t) by the formula 

I D (t) = U D (t)I D (0)U D \t). (51) 

The properties of the operator Io{t) are very similar to the properties of corresponding integrals of 
the motion in the case of Schrodinger equation and, of course, these Dirac invariants-operators tend 
to the integrals of the motion of the Schrodinger particle in nonrelativistic limit. Some theorems 
about Schrodinger integrals of the motion may be transfered to the case of the Dirac equation. 
For example, a function of the Dirac invariant-operator is another Dirac invariant-operator. This 
statement means, that if we have a function /(/^(t) ), this function may be always represented in 
the form 

f(I D (t)) = U D (t)f(I D (0))U D l (t). (52) 

But this formula follows from the obvious property 

/ (u D (t) I D (0) U D \t)) = U D (t) f(I D (0) ) U D \t) . (53) 

For example, for any power series with respect to the operator iu(£) it can be checked and extended 
for more complicated function. From this property, it follows that a free Dirac particle has also 
the following invariants: 

Pod, x od, (Podx od ), Pod L(t) . 

The eigenvalues of the Dirac invariants-operators do not depend on time. It means that if we 
search for the solution to the equation 

I D (t)4>(t) = \(t)^(t), (54) 

where bispinor ip(t) is a solution to the Dirac equation, the function X(t) is constant and does not 



depend on time. This statement follows from simple checking. We have 

U D (t) I D (0) U D \t) U D (t) V(0) = X(t) U D (t) V(0) . (55) 

After simple algebra, it is obvious that \(t) is eigenvalue of the time-independent operator Id(0) 
and due to this it does not depend on time as well. 

The third statement which may be generalized from the Schrodinger example is that if we have 
a solution to the Dirac equation ip{t) and act on it by the Dirac invariant operator In{t), the new 
bispinor <p(t) such that 

<p(t) = I D (t)i{>(t) (56) 

turns out to be the solution to the Dirac equation, too. This statement may be proved by simple 
substitution. For the solution ip(t), we have the expression 

m = u D (t)m- (57) 

Then, we have 

I D (t) m = U D (t) 72,(0) U D \t) U D (t) V>(0) , (58) 



but the bispinor 



<p{t) = U D (t)<p(0), (59) 



where 

(p{0) = I D (0)i/>(0), (60) 

satisfies the Dirac equation due to the property of the Dirac propagator Unit). 

The matrix elements of the Dirac evolution operator Uz>(t) in any representation satisfy the 
relations analogous to (|36|) and (|37 ) for the free Dirac particle. In the coordinate representation, 



for any Dirac Hamiltonian if we found the operators xo(t) and po(t), these equations hold. 



4 Newton— Wigner Operator 

The integral of motion, which is an initial position operator, turns out to play important role 



in considering a nonrelativistic quantum system [^, |24| . Due to this for relativistic Dirac particle 



it is worthy to study an analogous integral of motion, which plays the role of initial position of 
the particle. For the Dirac particle, the position operator itself has some pecularities and it differs 
from the position operator of nonrelativistic particle. The most appropriate operator playing the 



role of the position operator of the Dirac particle was considered by Newton and Wigner |18| , |19| . 
We will construct the time-dependent integral of motion, which has the physical meaning of 
the initial position operator. It means that this integral of motion at the initial moment of time 
must coincide with the Newton- Wigner position operator. 



The Newton- Wigner position operator [18[ |19j has the form (see, for example, |20 | 



where 



and the projector 



Qu = A + (P) (1 + 7o) J-%— (i £-) J-?-A + (p) , (61) 

V po + m \ dp k l V po + m 



2 2,2,2 

P =Pl+P2+ P3 



A+tf) = l( 1 + f*±^). (62) 



Also we used for the operator of coordinate in momentum representation the notation 

d 
Xi = i-—, i = l,2,3. (63) 

dp k 

For free Dirac particle, one can construct 3- vector-operator Qk (t) , which is integral of motion, 
using the general construction and invariant found in previous section, i.e., 

Qk(t) = U D (t)Q k U D 1 (t). (64) 



Since the operators 

U D (t) Po U D 1 (t)=p (65) 

and 

U D (£) A + (p) Uf (t) = A + (p) (66) 

are integrals of motion for free particle, the integral of motion, which for initial time t = coincides 
with the Newton- Wigner position operator is 



Q k (t) = A + (p) (1 + 7o (t)) JQ±^x h (t) ,/^±^A+ (p) . (67) 

f Po f Po 

Here the integral of motion x fc (t) = x Q Dk(t) is given in coordinate representation by (^Sj)-(^) and 



in momentum representation by the same formulas Q5q)--(fEq) , in which the operator x is replaced 
by i d/dp . The integral of motion 



lo(t)=U D (t) lo U D 1 (t) (68) 



in explicit form is 
7o (0 



cosh 2 £:r — £ 2 (m 2 — p 2 ) sinh 2 er — 2e 1 ap sinh err (e 1 sinh er + cosh er) 



i — 2e: 1 <rpsinh£:r (e x sinh er — cosh er) e 2 (?n 2 — p 2 ) sinh er — cosh er j 

(69) 

Thus, we get 3-vector-operator Qk if) , which is the integral of motion of free Dirac particle. This 

integral of motion "remembers" its initial maximally localized position in space. 



5 Short Pulses for Dirac Partilce 

Let us consider the model, in which the interaction term for Dirac particle is the term V (x, t) 



in fl4D|), i.e., we have 

H D = ap + f3m + V (x, t) . (70) 

We study what influence produces very short pulse of the potential. It is motivated by discus- 
sion [^] of heavy ion collisions so fast passing each the other in process of high energy reaction 
that the perturbation of the potential field lasts very short time in comparison with the periods 
of all motions in the system. In this case, the model of the potential V (x, t) may be described by 
^-function of time 

V (x, t) = kV (x) S (t - t') . (71) 

We study explicitly how the propagator for the free Dirac particle is changed if the potential 
of the type ([n]) switches in. This problem is interesting not only in the context of studying 
pair production in the process of charge collisions but also as the simple exactly solvable model 
of relativistic particle with nonstationary Dirac Hamiltonian. Then solving Eq. ( f48|) with the 
Hamiltonian (|7Dj), ( |7T| ) we can find that the propagator of the particle with S- type interaction in 
time may be factorized 

U D (if , t in ) = Ui (£ f , exp {-ikV) U f (£', t in ) . (72) 

Here U{ (£2, £1) is the propagator of free particle and the influence of <5-pulse at the moment t' 
gives the factor exp (—ikV) . The free propagator has been calculated in previous sections and is 
described by formula (|35|) . If we introduce the notation 

- i (t - *') = n ; -i (V - t ) = r 2 ; £1 = \fp\ + m 2 ; e 2 = \fpl + m 2 (73) 

and 

(exp [-inV]) pip2 = I exp [i (p 1 - p 2 ) x - inV (x)] dx , (74) 



in momentum representation, the propagator of the particle G (p\, P2, t, t ) with 5-pulse takes the 
form 



G (pi, p 2 , t, t ) = (exp [-inV]) J[ 



i=i 



( i i \ 

cosh £;Ti + £j m sinhejTj e^ <7 £>j sinh £jTj 

, e, rl a Pi sinhejTj coshe^ — e^ x m sinhejTj 

(75) 

That is exactly the form of the Green function for Dirac particle with 8 ( t — t J -time dependence 

in the potential, where t is the moment of 5-kick. For k — , the first term gives the 5 (j)\ — p^)- 

dependence and the product of two 4x4-matreces in ( [75] ) gives one matrix of the form fl35]). 

Let us find out the change of the form of the integral of motion for free Dirac particle related 

to the 5-kick, for which t' = 0. Obviously, for negative times the form of the integrals of motion 

operators is the same as for free particle. But for positive times the form is changed. So, for the 

form of the invariant, which is initial momentum, one obtaines 

(Po (*))«,' = 6 {-t)pS (p -p)+e (t)Jp" (exp [z«V]) w » (exp [iKV])* p y dp", (76) 

where 

6 (£) = 1 , t > ; 

9 (t) = , t < . 

The matrix (exp [— wcV]) p „ is proportional the unity matrix in four-space. 

Thus, we have found the kernel of the time-dependent integral of motion for the kicked Dirac 
particle. This integral of motion plays the role of the initial momentum of the relativistic quantum 
motion. 



6 Conclusion 

We have found some time-dependent integrals of motion for relativistic quantum particles. This 
consideration of time- dependent integrals of motion is motivated by the problem of pair creation 
in heavy ion collisions. Such processes may be modeled by fast time dependence of the interaction 
potential in Dirac equation. Due to this it is interesting to consider the constant of the motion in 
such type of models. The results of the paper are given in the explicit expressions (|38|)-(|42|) for the 
analog of the initial coordinate operator which is the integral of motion for the Dirac particle and 
in fl7"B] ) for the analog of the initial momentum, which is the integral of motion for the kicked Dirac 
particle. We have found the time-dependent integral of motion for free Dirac particle given by 
(67), which is initial position operator by Newton and Wigner. Also the time- dependent integrals 



of motion related to initial spin and angular momentum observables (|47|) of free Dirac particle 
are the result of the analysis presented. It is worthy to study other time-dependent integrals 
of motion, which may be related to the high-energy processes modeled by interaction potentials 
dependent on time. One should note that in all cases of the relativistic quantum problems, for 
which the propagators are known explicitly, it is possible to find the explicit expressions for the 
kernels of time-dependent integrals of motion using the procedure analogous to the procedure 
applied for finding the initial momentum invariant of the kicked Dirac particle. 
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